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The escape fraction at infinity is evaluated for massless particles produced in collisions of weakly 
interacting particles accreted into a density spike near the particle horizon of an extremal Kerr black 
hole, for the case of equatorial orbits. We compare with the Schwarzschild case, and argue that 
in the case of extremal black holes, redshifted signatures can be produced that could potentially 
explore the physics of particle collisions at centre of mass energies that extend beyond those of any 
feasible terrestrial accelerator. 
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I. INTRODUCTION 

Some of us have recently argued that rotating black 
holes surrounded by relic cold Dark Matter (DM) den- 
sity spikes may act as particle accelerators [1]. In the 
limit of maximally rotating, extremal Kerr black holes, 
we showed that collisions between particles, e.g. weakly 
interacting massive DM particles, may reach arbitrarily 
high center-of-mass energies. For some related works see 
[2~8]. 

In this sequel, we compute the escape fraction and flux 
at infinity of the highly blue-shifted particles. We demon- 
strate that the escape fraction is finite and of interest 
for any value of a, although the price one pays for sub- 
cxtrcmality is that the achievable Centre-of-Mass Frame 
(CMF) collision energy is limited. 

Our earlier work attracted considerable comment with 
regard to the fact that extremality for a Kerr black hole 
might not be achievable in nature [4, 5], although not 
all authors concur on this point [6]. Accretion spins up 
the hole but radiative back-reaction limits the degree of 
spin-up. The limit a < 0.9980 ± 0.0002 was derived [9] 
for a thin accretion disk, and criticised as being nonrc- 
alistic for more general situations [10]. More recently, 
the role of binary black hole mergers has been studied 
and superradiant extraction of angular momentum from 
the larger black hole has been shown to set the limit 
a < 0.9979 ± 0.0001 [11]. However we note that these 
limits apply only to astrophysical constraints, and that 
string theory may provide alternative options for gener- 
ating extremal black holes. 

Given these motivations, it is imperative that we ob- 
tain an estimate for the flux expected from the ultra- 
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energetic collisions that take place in the near-horizon 
region of Kerr black holes [1]. Concretely, we will fo- 
cus on the flux emitted by DM spikes that form around 
Intermediate-Mass Black Holes (IMBHs), as detailed in 
[12]. Such DM spikes are expected to extend down to 
distances that are of order the Schwarzschild horizon ra- 
dius, and to have a density that is essentially set by the 
annihilation rate of the DM. We take into account the 
fact that massless products of ultra-energetic collisions 
may be captured by the black hole by including a fully- 
relativistic escape function into the computation. 

The paper is organised as follows: we first set the scene 
in the next section, and present a relativistic formula for 
the total flux expected at infinity from collisions around 
an IMBH. This formula requires knowledge of the es- 
cape function, which encodes the proportion of massless 
products that arrive at infinity, as well as the DM den- 
sity profile around the black hole. We discuss the escape 
function in Section III, and the density profile in Sec- 
tion IV. We present our results in Section V, and end by 
expanding on some open issues that must be tackled to 
improve and extend the results of this paper. 



II. THE ANNIHILATION RATE, DENSITY 
PROFILE AND ESCAPE FRACTION. 

Our aim is to obtain the flux emanating from particle 
annihilations in a DM spike collapsing around a rotating 
black hole. We first explain how we intend to compute 
the flux, for steady-state (time-independent) conditions. 
The first ingredient is the DM number density in the 
vicinity of the black hole. We model this using a distribu- 
tion function in phase space n(x, v), where x and v are the 
coordinates of real space and velocity space respectively. 
This distribution function n(x, v), which is a Lorentz 
scalar, enters the annihilation rate quadratically. Defin- 
ing a(x, Vi,V2) as the annihilation cross section for a pair 

of DM particles and v re i(x, v\, v 2 ) = \j 1 — ( — u"u v ) 
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as the relative velocity of two DM particles with veloci- 
ties U\ and U 2 , we propose the following formula for the 



where dV\^ denotes the velocity space of the colliding 
particles. The preceding equation tells us the number 
of annihilations taking place locally. Our objective is of 
course to find the flux of massless particles escaping the 
gravitational field of the black hole. Only massless parti- 
cles produced with momenta satisfying certain conditions 
will propagate out to infinity. To obtain the number of 
escaping massless particles from the annihilation rate, we 
must therefore dress the rate of Eq. (1) with an "escape 
function" which encodes the probability that the massless 
annihilation products escape to infinity. We expect this 
function to go to zero at the horizon (such that all of the 
annihilation products get captured there) and approach 
one at infinity, so that all of the annihilation products 
very far from the black hole escape to infinity and con- 
tribute to the flux. Call the escape function e(x,Vi,v 2 ), 
where we make explicit its dependence on the velocities of 
the colliding DM particles. Folding this in with the anni- 
hilation rate above, we have that the number of massless 
particles escaping the black hole per unit time is 
given by 

iVoo = / / / e{x,v ll v 2 )dA^)d i x. (2) 
JVi Jv 2 Jv 

Note that in principle we can use a similar formulation to 
obtain the flux spectrum, but our sole purpose is to show 
that the flux is potentially observable. If the observer is 
a distance D away from the black hole, then to find the 
flux received per unit area per unit time by the observer 
we must divide the above by 4irD 2 as usual. Note that in 
this work we limit the phase space integration domains 
Vi,2 to in-falling DM particles moving on geodesies that 
are captured by the black hole, the so-called "plunge" 
orbits, as those are the ones which lead to collisions with 
high CMF energies. For our purposes, the spatial integral 
will be computed in a spherical volume V around the 
black hole, extending to a radius r/, measured in units 
where the Schwarzschild radius R s — ttibhG/c 2 is set 
to 1. The upper limit of the volume integral is a free 
variable (in principle), and we may choose to restrict it 
in order to focus on the size of the flux emerging from 
the near horizon region, where the CMF energies of the 
collisions are high. 

A crucial ingredient required to calculate the flux is the 
density profile of DM around the black hole. In this pa- 
per we use a density profile corresponding to the "spike" 
proposed in [12, 13]. This density profile includes the 
effects of DM annihilation on the DM density (as well 
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as other effects) and we shall make the assumption that 
it provides a good description in all regions down to the 
black hole horizon. We stress that this density profile is 
also non-relativistic, and so its convolution with a fully 
relativistic (and exact) escape function is a source of 
potential inconsistency but this method provides a rel- 
atively straight forward way to estimate the density of 
DM particles around a black hole. Finally, the types of 
density spikes described in [13] correspond to density pro- 
files around non-rotating black holes. The more energetic 
collisions occur for rotating black holes, and therefore in 
a more detailed analysis aimed at finding exact numbers 
a relativistic density profile is needed. Again we stress 
that we are only interested at this stage in an order of 
magnitude estimate for the size of the emergent flux so 
we leave such refinements to future investigations. 

Regardless of how one estimates the density profile, 
in order to calculate the emergent flux we need to know 
what fraction of collision products escape the gravita- 
tional field of the black hole. To make the estimate as 
straightforward and as clear as possible, we only consider 
the simplified case of the annihilation of two DM parti- 
cles into two massless particles (e.g. two neutrinos or two 
photons). The escaping fraction of collision products is 
an r-dependent quantity. For the case of a Schwarzschild 
black hole the escape fraction as a function of r is a known 
result and a useful treatment can be found in [14]. 

When considering the case of plunge orbits, the escape 
fraction becomes more complicated due to the fact that 
the collision products are boosted towards the black hole. 
In the next section, we describe the calculation of the 
escape fraction as a function of r for both Schwarzschild 
and Kerr black holes, including the relativistic effects of 
the boosts. 



III. ESCAPE FUNCTION 

In this section, we present an analytic treatment of 
the escape fraction for massless particles. These parti- 
cles are produced in the collisions of the DM particles 
and are consequently highly boosted towards the centre 
of the black hole. The details of this boost are dependent 
on the momentum of the DM particles. This of course 
is determined by the angular momentum of the two DM 
particles. It is a relatively straightforward exercise to 
construct the escape function for a massless particle pro- 
duced in the gravitational field of a Schwarzschild black 
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hole [14]. Here we want to do a similar calculation but 
with the massless particles boosted according to the kine- 
matics of the initial colliding DM particles. 

Let us first examine the motion of massless particles 
in the Locally- NonRotating- Frame (LNRF), which is the 
frame in which an observer who rotates with the local ge- 
ometry does their measurements. Working in the equa- 



torial plane, the LNRF momentum of a massless particle 
is given by 

pLNRF _ fpLNRF pLNRF pLNRF^j ^ 

where the components take the forms [15] 



P? BRF {r,a) 
P r LHRF (r,a,a,b) 
P^ HRF (r,a,b) 




^3 + a 2( 2 + r y ! + 2 + _4_II^ r 3 + 2(fl _ 6 )2 + r(fl2 _ 6 2) 

a r 3 + a 2 (2 + r)-2ab : 

r 3 + a 2 (2 + r) - 2ab 



(4) 



where we are working in units where the mass of the 
black hole and the Schwarzschild radius R s are set equal 
to 1. In this frame we can establish the conditions on the 
massless particle's 3-momentum as a function of r such 
that the it escapes the gravitational field of the black 
hole. This is most easily performed by analysing the 
turning points (i.e. where r — 0) of the massless particle 
as a function of r and b. 

The equation of motion for r on the equatorial plane 

is 

{b 2 -a 2 ) + 2(6 -a) 2 



r = 1 - 



(5) 



Setting this equal to zero and solving for b as a function 
of r we find the two solutions 



b(r) = 



-2a ± vVr 2 - 2r 3 



(6) 



We can now use these solutions to tell us whether a 
massless particle produced at some distance r, with par- 
ticular values of b and a (a = ±1 and determines whether 
the massless particle is moving to larger (+1) or smaller 
(— 1) radii) will escape the black hole. 

To do this for general a analytically is tricky, and we 
choose to investigate two values of a, namely a = 0, the 
Schwarzschild limit, and the extreme Kerr case of a = 1. 
Other values of a give results that are numerically not 
too different from the analytic results presented below 
and we leave the details of a numerical treatment for 
arbitrary a to future work [16]. Our goal here is to find 
an estimate for the emergent flux rather than to give a 
precise prediction. 

Starting with the Schwarzschild case, a — 0, the two 
solutions of Eq. 6 take the forms 

&rV)=±-i- (7) 

Vr-2 



We will for the remainder of the Schwarzschild discussion 
omit the superscript a = 0, to avoid clutter. By plotting 
these solutions in FIG. 1 we see that for the positive 
solution there is a minimum and for the negative solution 
there is a maximum. Both these extrema occur at r = 3 
and have values b + (3) = 3^/3 and 6_(3) = — 3\/3- We 
can now read off the conditions on b as a function of r 
that determine whether the massless particle will escape 
or be captured by the black hole. 
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FIG. 1: A plot of turning points for Schwarzschild black holes 
in terms of 6 as a function of r (in units where R s = 1). The 
upper (red) curve is 6" =0 (r) whereas the lower (blue) curve is 
b a -T°(r). 



The conditions for the massless particle to escape are 



for 



r < 



3 : a = +1 and - 3\/3 > b > 3\/3 , 



(8) 



which simply reflects the fact that a particle produced 
in the region r < 3, with a value of b in between the 
minimum and maximum, will move towards infinity un- 
hindered by the potential barrier. If the particle were 
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produced with b > 3y/3, then it will hit the potential 
barrier and bounce back into the black hole. Similar rea- 
soning then yields the following additional conditions for 
escape, which apply in the region r > 3: 

r > 3 : a = +1 and 6_ (r) > b > b+(r); (9) 
a = -land 3^3 > b > b + (r); (10) 
a = -1 and 6_(r) > 6 > -3^3. (11) 

For the cr = — 1 conditions, the massless particle will 
initially be moving to smaller r but will reach a turning 
point (at a distance r determined by the value of b(r)) 
and will be reflected by the effective potential leading to 
the escape of the massless particle. 

To calculate the escape function, we apply the follow- 
ing logic. We assume that the annihilation of two DM 
particles proceeds as XX ~~ ^ P1P2, where we assert that 
the piS are massless particles. As an estimate designed 
to provide a reasonable order of magnitude measure, we 
only concern ourselves with the question of whether pi 
escapes. By doing this, we avoid the complications intro- 
duced by trying to consistently analyse whether either, 
both or neither p\ or p 2 escape. 

We construct the measure by assuming that in the 
CMF of the collision, pi is produced isotropically, that 
is, there is no preferred direction. In the CMF, there is 
an r dependent range of angles that correspond to par- 
ticular values of b and a (evaluated in the LNRF) that 
satisfy the conditions for pi to escape. In other words, 
the conditions on the values of b and a for the massless 
particle to escape in the LNRF (Eqs. (8) - (11)) can be 
mapped onto a range of angles in the CMF of the colli- 
sion in which pi is produced. Because of the assumption 
that the production in the CMF is isotropic, the escape 
fraction can be computed in terms of the ratio of this 
range of angles over the full 2ir in the CMF on the equa- 
torial plane. The situation can therefore be summarised 
as follows: two DM particles collide at a given r with 
two respective angular momenta. We have two frames: 
the LNRF frame, which we can think of conventionally 
as the lab frame, and the CMF of the collision. The two 
frames are related by a boost, which is defined by the 
radial position and the angular momenta of the colliding 
DM particles. It is this boost that provides the mapping 
between the escape conditions on b and a and the range 
of angles in the CMF. 

One immediate important point is that by assuming 
the massless particles are only produced and escaping in 
the equatorial plane, we are potentially overestimating 
the escape fraction. However, this is not entirely clear as 
some particles could still escape even if they have trajec- 
tories that take them off the equatorial plane. For sim- 
plicity we restrict our analysis to the equatorial plane. 
In a paper to appear, a numerical analysis is performed, 
[16], where it is shown that this approximation is indeed 
a good one to the accuracy to which we wish to work in 
this analysis. 

To find the r-dependent range of angles in the CMF 



that corresponds to the full range of escape conditions 
for b and a in the LNRF, we can construct dot products 
between the 2 spatial components of the massless particle 
momentum vectors evaluated in the CMF (i.e. the r and 
9 components) corresponding to the extreme values of b 
for the appropriate values of a and solve for the angle 
between them. This angle is then the escape range in 
the CMF, and all massless particles produced inside it 
escape. 

First we need to find the expressions for the momenta 
in the CMF. To do this we simply apply a boost to the 
LNRF momenta for the massless particles such that 

pCMF _ ^ pLNRF 

is the CMF momenta and where A is a general Lorcntz 
boost in the equatorial plane that can be written as 

A = 

7 — /?7COSx — sinx/?7 

— cosx/?7 1 + cosx 2 (~l + 7) cosxsinx(-l + 7) 

— sinx/?7 cosxsinx(— 1 + 7) 1 + sinx 2 (— 1 + 7) 

where x, 7 and (3 are the usual parameters of the Lorentz 
boost. The exact form of this boost is determined by the 
momenta of the DM particles that are colliding together. 
Parameterising the momenta of the massive colliding par- 
ticles (where we have set the mass of the DM particle 
fi = 1) as 



Qi(r,ai,a) = {Qii,Qif,Qi§) 



(13) 



where the momentum vectors are parameterised by the 
angular momentum a, and we use Q instead of P to 
make it clear that this momentum vector is for a massive 
particle as opposed to a massless particle. The explicit 
forms after taking a = are [15] 

Q if (r,a = 0) = -pL=, (14) 



V 1 ^ 



Q i ?(r,a i ,a = 0) = _V ^ + 2a^ra ?) > ^ 



Qi§(r, a 4 ,a = 0) 



(16) 



The boost parameters in terms of these momentum com- 
ponents are then found to be 



/3 - 



((Qif + Qjf ) 2 + {Q ie + Q 3 e)'- 

(Qa + Qjt) 
Qa + Qa 



1/2 



cosx 



sinx 



{Q i{ + Q^) 2 - {Qir + Q 3 r) 2 - (Q i§ + Q 3§ ) 2 

Qir Qjf 



A2 



{Qir + Qjf) 2 + {Q i§ + Q jS ) 
{Qif + Q j f) 2 + {Q i§ + Q, § ) 2 ' 
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We are now in a position to construct the dot products. 
Using the boosted momenta of the massless final state 
particles we can construct the dot products. Remember 
that this is in the CMF now. Explicitly we construct the 
vectors as 



P(r, b, a) = P r CM (r, b,a,a = 0), P^ M (r, b, a, a = 0) 



With reference to the conditions set out in Eq. (8) for 
r < 3 the range of escape angles, A(j>, in the CMF reads 



Escape Fraction EF 
0.8 r 



A0 r<3 (r) 



arccos 



P(r 


3^3, 


+i; 


•P(r, 


-3^3, +1) 


P(r, 


3\/3, - 






P(r, 


-3^3, +1) 



with the escape fraction for r < 3, EF r<3 , given by 

A0(r) r<3 



EF r<6 {r) 



2tt 



(17) 



P(r, b+(r), +1) 


.P(r, 6_(r), +1) 


P(r, b+(r), +1) 


P(r, 6_(r), +1) 



+ 



arccos 



arccos 



P(r 


Mr), 


-1) 


.P(r, 


3V3, - 


-1) 


P(r, 


6+(r), 


-1) 


P(r 


3^3, - 


-1) 


P(r, 


-3\/3, 


-1) 


■P(r, 


6_(r), 


-1) 


P(r, 


-3V3, 


-1) 


P(r, 


6_(r), 


-1) 



with the escape fraction for r > 3, EF r>3 , given by 

Af>(r) r>3 



EF {r) = 



2tt 



(18) 



Putting this all together, the total escape fraction is 

^[A^< 3 (r)6(r~3) 
+A#>j»(r)(l-e(r-3))], (19) 



EF a=0 (r) 



We again state that we are only considering motion in 
the equatorial plane. The corresponding expressions for 
r > 3 are (following the conditions stated in Eqs. (9, 10, 
11)) 

A0 r>3 (r) = 



arccos 



where 0(r — 3) = 1 for r < 3 and zero otherwise. 

In addition to the dependence of EF on r there is an 
implicit dependence on the initial angular momenta of 
the two colliding DM particles a% and a.%. We can plot 
this escape fraction for any combinations of a.\ and ui 
and in FIG. 2 we show EF for a,\ = 2 and a 2 = —2 
as an example. With different combinations of angular 
momenta the size of the escape fraction does not change 
significantly. 

As a consistency check, we can remove the boost from 
the calculation and consider massless particles produced 




FIG. 2: Escape fraction for a massless particle in the gravita- 
tional field of a Schwarzschild black hole as a function of r for 
DM particles colliding with angular momentum ot\ = 2 and 
a 2 = -2. 



at some position r and allowed to have initial momenta 
in any direction in the LNRF. This is obviously a much 
simpler problem. Comparing with the calculation of the 
escape cone outlined in [14] we find complete agreement. 

We now move to the extremal Kerr case. The calcu- 
lation for Kerr follows the same steps outlined for the 
Schwarzschild case. Let us repeat the reasoning for clar- 
ity: two DM particles collide at a given r with two respec- 
tive angular momenta. We have two frames: the LNRF 
frame, which we can think of conventionally as the lab 
frame, and the CMF of the collision. The two frames 
are related by a boost, which is defined by the radial 
position and the angular momenta of the colliding DM 
particles. Isotropic production of the massless particle in 
the CMF is again assumed. The conditions for escape, 
given by a range in b and a, as observed in the LNRF, are 
once again computed by considering the turning points 
of the massless particle. The boost is then applied to 
convert the momentum of the escaping particle to the 
CMF. This allows us to map the escape conditions in the 
LNRF, which are conditions on b and <j, to a range of 
angles in the CMF. This is then divided by 2ir to obtain 
the escaping proportion of massless particles. Although 
the formula for the turning points is now more compli- 
cated, the reasoning is identical to the Schwarzschild case 
above. 

The first step is to analyse the turning points of r once 
again. Using Eq. (6) and setting a — 1 we have the two 
solutions 



bl =1 (r) = r- 
b a _ =1 (r) = -'- 



1. 



(20) 
(21) 



We have plotted these two solutions in FIG. 3. Note that 
the conditions for escape are now more complicated, but 
the logic is the same. If the particle has a = +1 and it 
does not encounter a potential barrier, then it escapes. 
If it has a — — 1 and it encounters a potential barrier, 
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FIG. 3: Plot of turning points for an extremal Kerr black 
hole in terms of b as a function of r. The dashed (red) line is 
&+ =1 (j") whereas the solid (blue) curve is 651 =1 (r). 



then it bounces off and escapes. The solution, 6" =1 (r), 
has a maximum at r — 4 with value 6° =1 (4) = —1. The 
solution 6 t J =1 (r) has a lowest value of 2 which occurs at 
the horizon, r — r^. Again, we will drop the superscript 
a = 1 from what follows, as for the remainder of this 
section we will be concerned only with extremal Kerr. 

From this plot, we can read off the conditions needed 
for the massless particle to escape (note that here 6° =1 (r) 
refers to the lower (negative) branch of the curve, as dis- 
played in FIG. 3). They are for 

r<4:cr = +l and - 7 < b < 6+(r); (22) 

a = -1 and 2 < b < b+(r), (23) 

and for 

r > 4 : a = +1 and b_(r) < b < b+(r); (24) 

a = -1 and 2 < b < b+(r); (25) 

a = -1 and b_(r) <b < -7 . (26) 

Let us give the physical meaning of the escape conditions 
for r < 4. In this particle which is moving away 

from the horizon (a = +1), will escape unless it hits the 
lower solid curve (6_(r)) of FIG. 3. To avoid that, it must 
have b > —7, which is the maximum of that curve. On 
the other hand, if it has a value of b that is larger than 
the dashed line (6+(r)), then it will hit that line and 
bounce back to the horizon: this explains the condition 
with a = +1. For r < 4 and a = — 1, the particle which 
is moving towards the horizon, then escape will occur if 
it hits the red line, thus having b+(r) > b > 2. Similar 
reasoning yields the remaining conditions. We now need 
our final state massless particle momenta in the CMF. We 
follow the same procedure as before but now the initial 
DM momenta have the form 



Qif(r,a = l,o,) = 
Qi§( r > a = l -> 0L i) = 



C -2ati 
frC{r - 1) 



- Aa t - (r - 2)c 



y/r{r - 1) 



a, 



r 

C' 



(27) 



where C — r 3 + r + 2. 

The boost parameters are now set by the choices of 
a of the initial DM particles. Using the components in 
Eq. (27) and the boost, we can calculate the CMF mo- 
menta of the massless final state particles. These CMF 
momenta will be functions of b, r and a and we can apply 
the conditions necessary for escape listed in Eqs. (22, 23, 
24, 25, 26). 

Following again the procedure outline above we now 
want to construct dot products of the vectors 



P^(r,b,a)= Pf M (r,b,a,a = l), P% M (r, b, a, a = 1) 



evaluated with parameters corresponding to the extreme 
escape values quoted in the escape conditions. The prime 
notation is there to remind the reader that we are evalu- 
ating this with a = 1 (not to be confused with the tilde 
notation for the a = Schwarzschild case). Thus, we 
are able to obtain the angular range in the CMF where 
produced massless particles escape. This is in complete 
analogy to the Schwarzschild case above. 

The range of escape angles that corresponds to the 
range in the values of b and a for r < 4, as listed in 
Eqs. (22, 23,) A^j , in the CMF reads 



A#S(r) = 



' P'(r, 6+(r), +l)~P'(r, ~7, +1) 
_\P'(r, MO, +l)\\P'(r, -7, +1)| 
P'(r, b+(r), -l).P'(r, 2, -1) 
\P'(r, 6+(r), -l)\\P'(r, 2, -1)| 



(28) 



For r > 4, applying the conditions in Eqs. (24, 25, 26) 
A0a=i takes the form 



AC5(r) 



P'(r, b+(r), +l).P'(r, 6_(r), +1) 



\P'(r, b+(r), +l)\\P'(r, 6_(r), +1)| 
" P'(r, Mr), -l)-P'(r, 2, -1) 
.|P'(r, b+(r), -l)\\P'(r, 2, -1 
P'(r, -7, -l).P'(r, MO, -1) 
\P'(r, -7, -1)| |P'(r, 6_(r), -1) 

The total r dependent escape function is given by 
1 



(29) 



EF a=1 (r) 



2tt 



[AC<l(r)6(r-4) 



+A^>t(r)(l-6(r-4))], (30) 
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FIG. 4: Escape fraction for a massless particle in the gravi- 
tational field of an extremal black hole as a function of r for 
DM particles colliding with angular momentum ot\ = 2 and 
a 2 = -2. 



where 0(r — 4) = 1 for r < 4 and zero otherwise. Again 
EF a= i(r) has an implicit dependence on the initial DM 
angular momentum a a and a^. In FIG. 4 we show a plot 
of the escape function as a function of r for initial angular 
momentum a a — 2 and a>b = —2. 



Using this escape fraction, we can not only calculate 
the emergent flux of massless particles but we can plot 
the escape fraction vs the CMF energy of the collision. 
Doing so gives a clear measure of how likely it is that a 
massless particle that is produced in a very high energy 
collision will escape the black hole. We are able to do 
this because the CMF energy of the collision is purely a 
function of the radius and the angular momenta of the 
colliding DM particles, [1], 



(<T) = r(r 2_2r + a 2) ( 2fl2(1 + r) _ 2a( " 2 + ai) - a 2"i(-2 + r) + 2(-l + r)r 2 



- J2(a - a 2 ) 2 - a\r + 2r 2 J2(a - a ± ) 2 - a\r + 2r 2 . (31) 



Given this relationship, we can construct the plot shown 
in FIG. 5, which displays the escape fraction versus CMF 
energy for two DM particles colliding with angular mo- 
mentum ax = 2 and a 2 = — 2 and is evaluated in the 
extreme Kerr limit, a = 1. It is noteworthy that even 
for CMF energies that are close to 1000 times the DM 
mass, the escape function is almost 0.001, which is non- 
negligible, especially if these collisions are taking place 
frequently which of course depends on the DM density 
in this region as well as the size of the cross section 
for DM collisions. To get an estimate for the size of 
the escape fraction for higher energies, a good approxi- 
mation for the parametric behaviour plotted in FIG. 5 
is EF pa OJE^fi. This approximation is good for 
E cm > 10/i, any lower and this approximate paramet- 
ric dependence breaks down. 



IV. THE DARK MATTER DENSITY PROFILE 

Now that we have the escape fraction for a rotating 
black hole, we must turn to the DM density profile in the 
area immediately around the black hole, a region where 
reliable modelling of the DM dynamics is difficult. In 
what follows, we shall make some simplifying assump- 
tions that enable us to make a reasonable estimate for 
that density. We consider an IMBH, taken to be of ref- 
erence mass 1O 5 M 5 M0, as predicted to exist in the inner 
galaxy as failed building blocks for the central supermas- 



Escape Fraction (EF) 
0.010 r 



0.006 



0.002 



0.000 1 



Ecm/fi 



FIG. 5: The escape fraction versus CMF energy for two DM 
particles colliding with angular momentum ai = 2 and ct\ = 
—2 in the gravitational field of a Kerr black hole with a = 1. 
The energy is displayed in units of the DM mass fi. 



sive black hole, and possibly existing throughout the halo 
as failed satellites [12]. Here M is the mass of the Sun 
and M 5 is a dimensionless number parameterising the 
absolute mass of the black hole. We caution that this is 
a hypothesis, and there may be other routes to building 
supermassive black holes. 

In addition, we assume that the DM consists of mas- 
sive weakly interacting massive particles that froze out 
in the early universe. Doing so gives us rough ball park 
figure for the size of the annihilation cross section. Such 
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particles are cold and form a density spike via adiabatic 
contraction during the formation phase of the IMBH by 
baryonic dissipation. For our purposes the details of the 
spike or any other part of the DM profile are irrelevant 
if they concern distances bigger than 0(1) x R s . All we 
need take from the profile is that the rising DM density 
as one approaches the horizon is tempered by its anni- 
hilation, so that in the immediate vicinity of the black 
hole an "annihilation plateau" is reached. The (larger) 
radius at which the plateau is reached is unimportant for 
our calculation because it is much larger than O(10)R s . 
On the other hand, the radius up to which the plateau 
continues (the inner radius of the plateau), we take to 
be the black hole horizon. Therefore, we assume that 
the DM density is essentially flat from the black hole 
horizon all the way up to values of the radius which are 
O(10) x R s , where R s = GM B h/c 2 . We emphasize that 
the details of the DM density profile at such radii are 
irrelevant for us, as we are interested in the most ener- 
getic particles resulting from annihilations taking place 
at radii that are 0(1) x R s . The density in the annihi- 
lation plateau is given by p p i — m x /(av)t, where (<jv) is 
the thermal-averaged annihilation rate per unit volume, 
m x is the DM mass, and t is the time-scale for spike for- 
mation. For details of the annihilation plateau, we refer 
the reader to [12, 13, 17, 18]. 

Although we now have an estimate for the total DM 
density in the near-horizon region, we are still short of 
information as to the composition of the DM particles 
with regard to their angular momenta. This is of crucial 
importance because the energies of the emitted photons 
are set by the angular momenta of the colliding DM par- 
ticles, and the very energetic collisions can occur for DM 
particles in plunge orbits. So what we must do is obtain 
an estimate of the proportion of particles in the plateau 
which are in fact within the plunge range of the black 
hole. The complexity of the problem means that this is 
beyond the scope of the paper, and, given that we are 
only interested in an order-of-magnitude estimate of the 
flux emitted in the immediate vicinity of the horizon, we 
simply assume a flat distribution in angular momentum. 
In addition, we assume that for the particles colliding 
near the horizon, the proportions of plunge and bound 
geodesies are comparable and so we take the absolute 
value of the density of both to be of order the height of 
the plateau. 

We remind the reader that the correct angular momen- 
tum decomposition of the DM particles in the plateau 
is likely to be very complicated, especially in the near- 
horizon region. Our aim in this analysis is to estimate the 
size of the total emergent flux originating from near hori- 
zon collisions. It turns out that we are almost insensitive 
to the details of the angular momentum distribution as 
we are effectively integrating over the full range of an- 
gular momentum states. If we wanted to calculate the 
energy spectrum of the emerging massless particles the 
choice of the angular momentum distribution becomes 
more important with the potential to change the shape 



of the spectrum. 

V. RESULTS 
The flux emitted in the near-horizon region 

In this section we put together the ingredients com- 
puted so far: the density in the near-horizon region and 
the escape fraction to calculate the flux. Going back 
to Eqs. (1) and (2), for plunge orbits we can exchange 
the dependence on the velocities for a dependence on 
the angular momentum of the colliding particles. The 
phase space distributions, n, can in principle be recast 
as functions of x and a. The a dependence encodes the 
distribution of angular momentum states which in gen- 
eral changes with x. This x dependent distribution is 
unknown (as discussed above) and we make the approx- 
imation that we can factorize the a dependence from n 
(such that this distribution is constant in x). We now 
assume that the remaining distribution function in a is 
flat across the allowed range in a. In addition, the es- 
cape function e(r, 0,1,0.2) changes slowly with ai and so 
we can treat it as a constant with respect to the integral 
over the angular momenta of the colliding particles. 

Further to this assumption we remind the reader that 
the escape fraction has been calculated for the equatorial 
plane only. To do the full volume integral we need to 
know how the escape function changes as we go off the 
equatorial plane. Instead of changing the form of the 
escape function to account for this variation (which is 
a task most easily completed numerically) we assume a 
spherical symmetry such that we can use the same escape 
fraction for the full volume. 

We can convert the number density n(x) to a mass 
density p(x) by dividing by the rest mass m x . For our 
purposes we set this density equal to the plateau density 
p p i. The final form for the flux arriving at some distance 
D away from the black hole simplifies to a spherically 
symmetric volume integral (a consequence of the assumed 
spherical symmetry of the density and escape function in 
the near-horizon region) and can be written as 

it?;/? 3 f Tl 

*°>^L *M*.«*.<*W (») 

where av is to be taken as some typical annihilation cross 
section, m x is the mass of the DM particle, and D is 
the distance to the black hole. A further complication 
which we do not address here is that the cross section 
will in general be energy and therefore r dependent. As 
we wish to make a simple estimate of the flux we take 
the cross section to be constant in energy. A further 
simplification is that we assume that av is dominated by 
an s-wave term, which is velocity independent. In doing 
so av is assumed to be just a constant. In general there 
will be velocity dependent terms and consequently we 
would replace av with the equivalent thermally averaged 
cross section. 
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Note also that the form of the flux in Eq. 2 is similar 
to the flux spectrum constructed in Ref. [12] with the ad- 
dition of the escape function. Ref. [12] deals with the flux 
spectrum d$/dE which requires an additional ingredi- 
ent dNj/dE, which is the number of secondary particles 
produced per annihilation. We are interested in the total 
flux of particles, rather than the flux spectrum (the latter 
is obviously very important, but its calculation is much 
harder in a relativistic setting). We can trivially convert 
the flux spectrum to total flux by assuming that photons 
reaching us will have energies bunched around the mass 
of the DM particle, so that it is safe for our purposes to 
write dN^/dE = S(E — m x ) and integrate trivially over 
energy to obtain the total flux in Eq. (32). 

The escape fraction e(r, a.\,ot2) was computed in Sec- 
tion III. It is relativistic, and so its convolution with 
a non-relativistic density is one of the major approxi- 
mations in this work. As the escape function is slowly 
varying with on we may choose a\,a 2 as we wish (if we 
knew the angular-momentum distribution of the plateau 
then we would have integrated over these two vari- 
ables weighted by the distribution function in angular- 
momentum space). We choose one of ol\ and a 2 to be the 
critical angular-momenta corresponding to the highest- 
energy collisions. This choice will not alter the order of 
magnitude of the results below. The density p p i(r) is sim- 
ply the moss density of DM particles that are on plunge 
orbits. As outlined above, we will use p p i(r) = again 
implicitly assuming that the proportion of plunge parti- 
cles is of the same size as non-plunge particles. 

Note that we are integrating from the horizon up to 
an upper radius n which is a few horizon lengths. The 
factor of i?g is necessary for converting the integral from 
R s = 1 units to SI units. Plugging everything in we are 
left with 

R« m 



l(n,a) = / r 2 e(r, a 1 ,a 2 )dr , 

Jr h 

(33) 

where we have isolated the dimcnsionlcss integral I(r\ , a) 
inside which is found all of the dependence on the rota- 
tion of the black hole. This dependence on the rotation 
parameter a, which enters into the limits of the inte- 
gral as well as the integrand, is obviously rather com- 
plicated. Computing this integral for extremal Kerr and 
for Schwarzschild, with n = + 4 for both (and with 
ai = 2, ct2 = —2 for Kerr and a\ — 4, ot 2 = —4 for 
Schwarzschild), we obtain Ik ~ 20 and Is ~ 40. For 
comparison, we note that if we put 1 in place of the es- 
cape fraction we obtain the values Ik ~ 40 and Is ~ 70. 

We can now give an estimate of the expected emergent 
flux. We use the following fiducial values for our various 
parameters: we set the mass of the black hole using M 5 
and set it equal to M = 40 in units of 10 5 solar masses. 
We use (av)o = 10 _28 cm 2 s _1 , set the distance between 



us and the black hole to Dq = 10 kpc and assume that 
the growth time-scale of the black hole is to = 10 10 years. 
The total flux is then given by: 

(34) 

where <f>o = 3.41 km _2 ycar _1 . Note that by writing the 
flux in this fashion we have placed all of the dependence 
on the geometry into the integral X(ri,a). We see that 
IMBHs can be bright DM annihilation sources, and this 
has been explored and constrained in the literature [19]. 
For the case of extremal Kerr, integrating up to n = 5, 
we obtain the flux to be 70 km _2 year _1 . This is well 
within the reach of planned neutrino detectors. If we in- 
tegrate only up to n = 1.1 for extremal Kerr, thereby 
including the emergent flux of massless particles coming 
from the most energetic annihilations, we obtain an ex- 
pected number of 0.26 km _2 year _1 . Futuristic large area 
high energy neutrino/particle detectors in space such as 
OWL (http : I / owl.gsfc.nasa.gov/) would have an ef- 
fective detection area solid angle efficiency product in 
excess of 10 5 km 2 sr and could detect such signals over a 
few days (admittedly only above 10 19 eV). We note be- 
low that Penrose boosting may allow such energies to be 
achievable at infinity. 

While the actual flux from an IMBH depends on its 
distance from us and on the other reference parameters, 
which are uncertain by several orders of magnitude, we 
note that our calculation of a finite escape fraction which 
is weakly dependent on the black hole angular momen- 
tum means that even detection of a single event might 
open a window on new physics, attainable only via par- 
ticle collisions at extreme CMF energies. Our final result 
should be viewed as an order of magnitude estimate of 
the flux received on earth, which we find to be the same 
with or without black hole rotation. On the other hand, 
given that very high (even arbitrarily high) energies can 
be achieved in the annihilations around a rotating black 
hole, while the Schwarzschild geometry only allows max- 
imal CMF energies of 2^/5m dm , our results imply that 
the spectrum received from a rotating black hole should 
typically contain signatures of highly energetic products. 
A next step is to compute the spectrum as a function of 
a to confirm this expectation. We leave this for future 
work. 



VI. OPEN ISSUES 

We obtained above an order-of-magnitude estimate for 
the blue-shifted flux emitted by DM annihilations in the 
near-horizon region of an extremal Kerr black hole. We 
made several approximations and assumptions along the 
way that made the calculation tractable. We assumed 
that the escape function calculated purely on the equa- 
torial plane was sufficient, restricting the colliding DM 
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particles and the escaping particles to motion only on the 
equatorial plane. We also assumed that the DM "spike" 
around the black hole was non-relativistic and spherically 
symmetric, and that it was flat right down to the black 
hole horizon. In our view, the former of these two con- 
cerns is unlikely to be a source of large errors, because the 
escape fraction goes from zero to one regardless of the di- 
mensionality or the geometry. It is very easy to check us- 
ing our results above that the change from Schwarzschild 
to Kerr makes very little difference to the escape func- 
tion. Moreover, some work we have done, and upcoming 
work [16], shows that going to three-dimensions does not 
change the results significantly, thus corroborating our 
intuition. The concern about the DM density profile is 
rather more serious, because the density enters the flux 
quadratically. Several effects need to be studied more rig- 
orously, chief among them being the effect of the rotating 
geometry on the DM spike, as well as the feasibility of a 
more accurate modelling of the spike close to the black 
hole horizon in a fully relativistic fashion. The advan- 
tage of such a construction is that it should naturally 
yield the composition of angular momentum in the DM 
spike, which is crucial for the flux and its spectrum. 

We should mention here that a reliable computation 
of the flux spectrum requires another ingredient: an ex- 
plicit particle physics model of the DM and its interac- 
tions. This is necessary because it determines the energy 
dependence of the differential cross section of the an- 
nihilation diagrams, which, coupled with the geometry- 
dependent CMF energy of the collisions, feeds directly 
into the differential flux spectrum. The flux spectrum is 
therefore expected to be a far more revealing signature 
than the total flux, as it should contain clues as to the 
nature and composition of the DM spike. For example, 



we naively expect the spectrum to bunch around the DM 
mass, which in itself would be an exciting observation. 

There are further necessary refinements of our work 
that we would like to comment on. The computation of 
the CMF energy in the Kerr field reference [1] neglected 
back- reaction effects, which were studied in [4, 5]. This 
issue deserves further attention. Finally, supermassive 
black holes, the nearest one of which is at our Galactic 
Centre, should not be ruled out as possible sources 
of annihilating DM. Even after partial destruction by 
a complex merging history, something which is not 
necessarily an issue very close to R s , cusps, albeit with 
softer profiles , are regenerated by dynamical processes 
[20]. Finally, one of the most intriguing aspects, that 
we reserve for future work, is the exploration of the 
Penrose process [21], which may allow the escaping 
annihilation products to tap the rotational energy of the 
Kerr black hole, provided that the collisions occur within 
the black hole ergosphere. This effect was briefly consid- 
ered in [6], but the potential remains to be fully explored. 
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